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Abstract 

We propose the second moment of the Husimi distribution as a measure of complexity of quantum 
states. The inverse of this quantity represents the effective volume in phase space occupied by 
the Husimi distribution, and has a good correspondence with chaoticity of classical system. Its 
properties are similar to the classical entropy proposed by Wehrl, but it is much easier to calculate 
numerically. We calculate this quantity in the quartic oscillator model, and show that it works 
well as a measure of chaoticity of quantum states. 

PACS numbers: 05.45.Mt, 05.45.Ac 
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I. INTRODUCTION 



Although the quantum manifestation of chaos has been extensively studied over the past 
few decades, to define "quantum chaos" is still the main problem in this field. The direct 
extension of the definition applicable to classical chaos seems to fail because of the hnearity of 
the Schrodinger equation. Since quantum mechanics contains classical mechanics as a limit, 
however, there must be something in quantum mechanics that produces classical chaos in 
the classical limit. 

There have been many attempts to define a measure of quantum chaos. There are some 
measures using the level statistics, but much more information could be obtained from 
analysises of individual quantum states. We can classify measures of complexity of quantum 
states into two types: 

1. complexity of pure states 

2. complexity of an ensemble of quantum states 

For example, the von Neumann entropy of the density matrix belongs to type (2). In this 
paper, we focus on complexity of type (1). 

Some quantities defined in terms of the expansion coefficients are often used in numerical 
calculations. For instance, suppose a quantum state is expanded in an appropriate basis 

k> = Ecd^)- (1) 

i 

Let us define Pi — |cj|^. Then the information (Shannon) entropy 

>5 = -^p,lnp,, (2) 

i 

and moments of the distribution 

Mk=Y:pt (3) 

i 

are measures of localization (or delocalization) with respect to this basis. In particular, the 
inverse of the second moment Mg"^ is called the number of principal components (NFC): 
it becomes unity when the state has only one component, while it becomes n when the 
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probability is equally distributed over n basis vectors. Such quantities are easy to calculate. 
However, an obvious defect is that the definition depends on the basis. 

Wehrl has proposed a good measure of complexity of quantum states based on the Husimi 
distribution function Ph{PjQ) I31[0' called it "classical entropy", 



while in his paper "quantum entropy" denotes the von Neumann entropy of a density matrix. 
Note that the classical entropy is applicable and takes various values also for pure states, 
while the quantum entropy always vanishes for them. Although Wehrl introduced the clas- 
sical entropy as an approximation to the quantum entropy, their values are not necessarily 
close even in the limit h 0. (See the discussion in [Q.) According to our classification, we 
discuss the classical entropy to describe complexity of type (1), i.e. of pure states. 

Chaoticity in classical mechanics can be characterized by the delocalization of orbits. 
In integrable systems, there are many constraints from symmetries, such that orbits are 
confined to low- dimensional tori. As the system becomes chaotic, the tori are destroyed 
and orbits can spread to higher dimensional space. In highly chaotic systems, orbits spread 
uniformly over the equi-energy surface. Such systems are called ergodic. Among the several 
conditions to characterize chaoticity of classical mechanics, the ergodicity is a rather weak 
one. For example, any one-dimensional time-independent Hamiltonian system is ergodic, 
but never chaotic. However, in physically natural situations of many-dimensional systems, 
the ergodicity works as a definition of classical chaos. 

We can expect a similar behavior in quantum mechanics. The Husimi function is a 
function on phase space, and takes only non-negative values while the Wigner function 
can be negative and is usually violently oscillating Hence the Husimi function can be 
regarded as a probability distribution in phase space, and its order of delocalization can be 
a measure of chaoticity of quantum states. In that sense, Wehrl's classical entropy seems to 
be a good candidate measure. 

However, Wehrl's entropy is defined as an integral over phase space, and is not easy 
to calculate numerically. To the best of the authors' knowledge, there is no calculation 
for a more than one-dimensional system. (A calculation for a one-dimensional time-driven 
system is in [^].) The difficulty mainly comes from the redundancy of the Husimi function: 
the Husimi function of k dimensional system is a function on 2k dimensional space, while a 




(4) 
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function on k dimensional space is enough to keep the same information in either coordinate 
or momentum representation. 

Complex pure quantum states are usually represented as sets of expansion coefficients 
with respect to a basis. This will be the case also in our following discussion. When we 
calculate Wehrl's entropy of a quantum state, the normal procedure is the following: 

1. Calculate the Husimi function pn on many sampling points {(Pj, Qi)}- 

2. Take an average (ph^t^Ph) over the sampling points. 

However, this procedure seems to be overly excessive because we have all the information 
about the quantum state in the set of expansion coefficients. There must exists some formula 
to calculate the average directly from the coefficients. The main concern of this paper is to 
avoid the redundancy of the Husimi function in numerical evaluation of a suitable measure 
of complexity. In other words, we wish to know how to get some averages related to the 
Husimi function without actually calculating the Husimi function itself. 

It seems difficult to derive a simple formula for the entropy because of the transcendental 
logarithmic function. However, a simple formula is possible for an algebraic function. The 
second moment, which is the average of the square of the distribution function, is especially 
easy to calculate. We will show later a formula in which the second moment of the Husimi 
distribution is expressed directly in terms of expansion coefficients in the harmonic oscillator 
basis. (If the quantum states are given by expansion coefficients in another basis, we should 
calculate the transformation matrix and change the basis to the harmonic oscillator basis.) 
The numerical effort is of order N"^ for a quantum state where N is the number of basis 
vectors. If we calculate the second moments for all quantum states given by diagonalization 
of a X matrix, the numerical effort is of order A^^, which is the same as the order of 
the diagonalization of the matrix. 

We represent the inverse of the second moment by W2: 



W2 represents the effective phase space volume occupied by the Husimi function. For exam- 
ple, if ph takes the same value over a region with volume V and takes zero value outside of 
it, W2 — V. The unit of W2 is the Planck cell volume. 



W2{ph) 



1 



(5) 



M2{pHy 




(6) 
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Next we summarize the main points of this paper. We use the second moment of the 
Husimi distribution to define complexity of quantum states. It is a measure of delocahzation 
of the Husimi distribution in phase space. It is defined in a base-independent way and has 
a good correspondence with complexity of classical mechanics. Moreover, we can calculate 
it directly from expansion coefficients, without calculating the redundant Husimi function. 
Therefore it is not so difficult to evaluate numerically. 

This paper is organized as follows: In section |T| we summarize the definition and some 
properties of the Husimi distribution function. In section |T| we derive the formula to 
calculate the second moment of the Husimi distribution directly from expansion coefficients 



in the harmonic oscillator basis. In section IV we introduce a model Hamiltonian and show 



numerical results which illustrate the meaning of W2 as a measure of complexity. The final 
section is devoted to a summary. In appendix we show details of semiclassical calculations 
of the second moment of the Husimi distribution in integrable and ergodic limits. 

II. HUSIMI DISTRIBUTION FUNCTION 

In this section, we review some properties of the Husimi function [||. We restrict our- 
selves to a one-dimensional system for simplicity, but the generalization to many- dimensional 
systems is straightforward. 

The Husimi function of a quantum state \(p) is defined as 

phAp,q) = \{^,M^)\'- (7) 

Here, \z, A) is a coherent state defined as an eigenstate for a complex eigenvalue z 

ax\z, A) = z\z, A), (8) 
where a\ is an operator with A as an arbitrary parameter 

The real and imaginary parts of z are related to the phase space point {p, q) by 

q=S^{z + z), (10) 



It is not difficult to derive the relation 



X exp 



(12) 



where pw is the Wigner function 

Pw{p.q) = I dr]{q-r^/2\^){^\q + r]/2)e'^^/\ (13) 

A state \q) represents an eigenstate of the coordinate representation. There we see that the 
Husimi function is a kind of the coarse-grained Wigner function. The main advantage of the 
Husimi function is that it is non-negative, as is obvious from the definition (|^). Hence the 
Husimi function can be formally regarded as a probability distribution, whereas the Wigner 
function can not. 



III. SECOND MOMENT OF THE HUSIMI DISTRIBUTION 



In this section, we consider the second moment of the Husimi distribution 

M.ipH) = I^PHiP,qr, (14) 
— \{z\^)\\ (15) 

Here and in the following we omit A for simplicity. 

Suppose the quantum state \ip) is represented in the harmonic oscillator basis 



= IIc„|n), (16) 



ra=0 



then M2 can be expanded as 



M2{Ph) = E E """"""^""^ / ^e-2|.p^n+n'^m+,n'^ (^7) 



n,n' m,m' Vn\n'\m\m'\ J TT 



n,Ti' m,m' 



Here, we used the following formulas 



(19) 



{n\z) = 6-1^1^/^4=, (20) 
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and 



—e-\'\z^z'^ = n\5nm- (21) 

TT 



As a result, we obtain the formula 

oo 



M2{ph) = \ Y.\Bl\\ (22) 



2 

^ L=0 



where 



L 



Since there is a truncation < in the expansion (plGf ) in any real calculation, the sum 
(P^) is also finite. The numerical effort is of order A^^ for a quantum state. 

The generalization to many-dimensional systems is straightforward. The result is 

M.ipH) = 2-^$:|s^f , (24) 
L 



where k is the dimension of the system and 



% = E 



c-cr (25) 



Here, L and j are fc-dimensional vectors whose components are non-negative integers. The 
factorial of those vectors means the product of all factorials of the vector components. For 
example. 



k 



L\ = Y[U. (26) 



i=l 

The absolute value denotes 

k 

\L\=T.L^, (27) 

1=1 

and j < L means ji < Li for 



IV. NUMERICAL RESULTS 
A. Model 

The model Hamiltonian we are considering here is that of a two-dimensional quartic 
oscillator 



H=\{vl+ VD + ^(^' + y') - kxV, (28) 



which has been adopted by many authors for the studies of level statistics and wave functions 
p. We put h = 1 and regard all quantities as dimensionless. 
This model has the simple scaling property, 

H{a^p^,a'^Py,ax,ay) = a^H{p^,py,x,y), (29) 

which means that the energy of the system does not essentially change the dynamics. How- 
ever, the parameter k changes the nature of the system. The system is separable and 
integrable at = 0, and becomes chaotic as k increases. It is unbound for k > 1. Meyer 
^ has shown that for large k values (> 0.4) the classical phase space structure is almost 
completely chaotic. 

This system has a discrete symmetry called C^^: The Hamiltonian is invariant with 
respect to reflections about x-axis, y-axis and also about the line x = y. In this paper we 
treat only quantum states which are symmetric under all reflections. (This symmetry class 
is labeled Ai in 0.) 



B. Structure of the classical phase space 

Fig.|l| shows the Poincare surfaces of section of this system at A; = 0.0, 0.2, 0.4, and 0.6. At 
k = 0.0, the phase space is completely occupied by tori. These tori are partially destroyed 
at A; = 0.2, and most of them seem to vanish at k = 0.4. However, there are still many 
islands of stability. The two most significant islands are at x = 0. These correspond to the 
stable linear orbits on x = y and x = —y. At k = 0.6, we can not recognize any structure. 
In this case, more than 90% of the orbits are unstable according to IQ. 



C. Diagonalization of the Hamiltonian 

We diagonalize the Hamiltonian (|28D in the harmonic oscillator basis. The procedure we 
used is essentially the same as that of Zimmermann et al. [§. The bases belonging to Ai 
can be written as 

K,n,)^, = ^M±hl^, (30) 
V2(1 + 5„,,„J 

where Ux and Uy are even non-negative integer and > Uy. Matrix elements of the 
Hamiltonian (pHf ) can be calculated analytically, and we obtain eigenenergies and eigenstates 

8 



by the diagonalization of this matrix. The truncated Hilbert space is spanned by the basis 
vectors with < rix + riy < 270, whose dimension is 4692. 

Among the 4692 eigenstates we obtain by the diagonahzation, those having very large 
energies are not rehable because of the truncation error. Since the range of the vahdity 
very much depends on the oscillator frequency u of the harmonic oscillator basis, we should 
choose uj to optimize the diagonalization. Because of the variational principle, eigenvalues 
obtained in the restricted Hilbert space are always higher than the real values, and the level 
density in that space is always smaller than the real one. Therefore the minimization of trH 
(the sum of the eigenvalues) can be a criterion to choose the optimal u. 

We put UJ = 7.0, which is chosen roughly to optimize the diagonalization for k = 0.6, for 
all k. The comparison between the obtained level density and the semiclassical one shows 
that the maximum reliable energy is E = 800 ~ 1000. 

D. The second moment of the Husimi distribution and the number of principal 
components 

Figs. D, ^, ^ and § show the results of the numerical calculation of the inverse of the 
second moment {W2) of the Husimi distribution for k = 0.0, 0.2, 0.4 and 0.6. We plotted 
them for energy eigenstates with E < 600. In each figure, 1000 ~ 1500 eigenstates are 
plotted. It takes about 30 minutes on the NEC SX-5 in Osaka University to calculate the 
second moments for all (4692) eigenstates. 

When we calculate the Husimi distribution, A in (0) is a free parameter. The simplest 
choice is to set A equal to the harmonic oscillator frequency of the basis used for the diag- 
onalization. In this case, we can directly calculate the second moment using the formulas 
(^) and (P^l). The results in Figs. 0, 0, ^ and ^ are calculated using the optimized value 
A = 7.0. 

We also calculated W2 for different values of A. Fig. || shows the results for A = 4.0, 7.0 
and 10.0 at A; = 0.2. Qualitative features of the figures seem unchanged, at least when 
A is not so far from the optimized value. When A is far from the optimized value, it is 
hard to obtain reliable results because the number of basis vectors we need to represent the 
eigenstates is huge. 

At A; = 0.0, the system is separable into two one-dimensional systems. Therefore eigen- 
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states of the original system can be specified by quantum numbers rrix and my, wliich label 
eigenstates of the one-dimensional systems. 

We can assume > rriy without loss of generality in the class Ai. Fig.|^ shows the 
results of a semiclassical calculation based on the torus quantization. (For details of the 
calculation, see Appendix [A 1| .) Semiclassical results and full quantum results have a good 
correspondence, except for the cases where ~ rriy and rriy ~ 0. When rrix ~ my, the 
assumption (|A21D does not seem to apply. When my ~ 0, the quantum number is too small 
to use the semiclassical approximation. 

In Fig.|[ the regular structure of Fig.|^ is partially destroyed, but there are still many 
regular series of eigenstates. The most significant series, which is in the lowest part of the 
Fig.^], corresponds to the stable diagonal periodic orbits ai x = y and x = —y. Fig.H shows 
an eigenstate in the regular series. There are also some series of regular eigenstates which 
are excited in the transverse direction to the periodic orbits. 

The lowest series in Fig.§ corresponds to the periodic orbits at x = and y = 0, which 
are stable at = 0.0. However, these orbits are unstable at k > 0, and this series seems to 
disappear as k increases. 

The eigenstates corresponding to the torus with = Ey are in the middle of Fig.^ (See 
also Fig.|^.) As k increases, this torus is quickly destroyed and two stable diagonal orbits 
are left. Islands of stability around these two orbits become the most significant structure 
in this system, and the eigenstates localized around these orbits form the lowest series in 
Fig.|. 

Aik = 0.4, as seen from Fig.|^, there are still regular series that correspond to the diagonal 
orbits and the first excited states in the transverse direction. However, most structures seem 
to have been destroyed and many eigenstates are near the ergodic limit, as based on the 
Berry- Voros hypothesis. (See Appendix |A 2| .) 

At A; = 0.6, no clear structure can be seen in Fig.||, and the values of W2 go up as a 
whole. However, the ergodic limit based on the Berry- Voros hypothesis is still not reached, 
and some eigenstates have much smaller values than the limit. 

There may be several reasons for this. One reason we can think of is that the Wigner 
function in the Berry- Voros hypothesis does not satisfy the pure state condition of the density 
matrix fp' = p. Therefore we might obtain better estimates of W2 by taking into account 
the pure state condition by successive iterations. (See chapter 8 in 0.) Another reason. 
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which seems more important, is that there are weak locahzation phenomena hke scars [jlO 
We found many scarred eigenstates with small W2 at k = 0.6. 

At low energies {E < 50), W2 of some eigenstates reach the ergodic limit irrespective of 
the value of k, but the figures of these wave functions does not seem chaotic. The reason 
for this is probably that the semiclassical approximation we used here is not good in this 
region. Our approximation is based on the idea that the local volume occupied by the 
Husimi distribution is the volume element of the invariant manifold (invariant torus or equi- 
energy surface) multiplied by the thickness factor of order h''^'^, where / is the codimension 
of the manifold. (See Appendix ^.) However, this approximation fails when the scale of the 
invariant manifold is comparable to the Planck constant. For example, in the low energy 
part of Fig. the ergodic limit line is under the limit based on the torus quantization. There 
seems no clear distinction between chaotic and regular eigenstates at very low energies. 

In Fig. 1^, we plotted the NPC at k = 0.0, 0.2, 0.4 and 0.6. We can see regular structure 
in Fig.^ (a). The structure is gradually broken as k increases, and the points are lifted as a 
whole. In that sense, the behavior of NPC is similar to that of W2- However, by comparing 
NPC and W2 of each eigenstate, we can see that the two quantities are not so similar. For 
instance, at k = 0.4, we can see regular series of eigenstates related to the diagonal orbits in 
both Fig. ^ and ^ (c). However, the values of the NPC are not so small, and there are many 
states lower than the series which seem to have no clear structure. Therefore, in general, 
W2 seems more reliable than the NPC as a measure of complexity of quantum states. Some 
eigenstates that have small NPC are related to the scars of the unstable linear orbits at 



a; = and y = 0. (See Fig. |T0[) Importance of these eigenstates in the response function 
has been reported in [|ll|]. 



V. SUMMARY 



In this paper we proposed the second moment of the Husimi distribution as a definition of 
complexity of quantum states. Its inverse (we represent it by W2) shows the effective volume 
occupied by the Husimi distribution function, and serves as a measure of delocalization in 
phase space. We calculated it for a quartic oscillator model, and showed that it has a good 
correspondence with chaoticity of the classical system. We can calculate it directly from 
expansion coefficients without numerical integration. Therefore the calculation is not so 
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time-consuming, and there is no numerical error except for that contained in the quantum 
state itself. 

In the integrable case {k = 0.0), the values of W2 have a regular structure. As k increases, 
the structure is gradually destroyed and the values go up near the ergodic limit line as a 
whole. At A; = 0.2 and 0.4, there are many regular series related to stable islands in classical 
phase space. Even at k = 0.6, the values of W2 of some eigenstates are much lower than the 
ergodic limit based on the Berry- Voros hypothesis, and they seem to be related to scars of 
unstable periodic orbits. 

To generalize the idea of this paper to many-body systems using generalized coherent 



states is interesting. It will be reported in the forthcoming paper |jT^. 
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APPENDIX A: SEMICLASSICAL CALCULATIONS 

In this appendix we give, for comparison, semiclassical estimates for the second moment 
M2. 



1. Integrable case (A; = 0) 

In this case, the model Hamiltonian is separable: 

HiPx,Py,x,y) = h{p^,x) + h{py,y), (Al) 

h{p,x) = ^{p' + x^). (A2) 

Therefore, we first derive a semiclassical formula for the Husimi function of eigenstates of 
(^). After that, it is easy to construct the semiclassical eigenstates of (|A1|) taking into 



account discrete symmetries. 
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a. One- dimensional problem 



The action variable of the one- dimensional system is 



27r 
4 

3^ 



pdx, 



(A3) 
(A4) 



C 



1 dx 



1 



(A5) 



/o y/l-X^ 4V27f V4. 
The semiclassical eigenvalues are determined by the Bohr- Sommerf eld quantization condition 



I{Em) = (2m + l)7ih, 



(A6) 



l/37r\^/3/ l\4/3 

Err, = - \ — ] I m + 



I.e. 



2 V4Cy V " ' 2 
The semiclassical Wigner function of the eigenstate is 

and the Husimi function is obtained by Gaussian smearing thereof: 



PH,x{p,q) 
1 



dp dq exp 



-1 Jw f^2 , iP-P' 



/\2 ■ 



Pw{p,q')- 



(A7) 



(A8) 



(A9) 
(AlO) 



By changing coordinates q q/\^,q' q'\/X,p — > \/Xp,p' — > -\/Ap', this becomes equiva- 
lent to using the modified Hamiltonian 



H=-p' + 

2^ 2A2^ ' 



(All) 



and putting parameter A = 1 in the Gaussian. 

In this case, the Husimi function near the equi-energy surface can be written approxi- 
mately PI 



2Vnh 



T\gradH{s) 



exp 



(e-^o 



(A12) 
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Here, s is a coordinate which parameterizes the equi-energy surface. ^ is the other coordinate, 
and ^ = ^0 corresponds to a point on the surface, s and are assumed to be orthonormal. 
T is the period, whose exphcit form is 

T = 4C{2E)-^/\ (A13) 

To calculate M2{ph), we introduce new coordinates {E,9) as 



Then 



2E 

p = ]l — sm9, (A14) 
q = {2Ef''VXcos^/^e. (A15) 



M2{Ph) = l'^PHAi.s)\ (A16) 



ds 



V2Tih 

2^2 



d9 



dd 



1 



\gradE{d)\'^' 



1 n-Kfi r/^ de 



32C2V E Jo ./icos4^ + -^sin2^cos0"' 

V 4\/2ii; 



(A17) 
(A18) 



b. Two-dimensional solutions 



Since the model (|28|) has some discrete symmetries, we have to take them into account 
when we construct semiclassical eigenstates. For example, in the class Ai, an eigenstate is 
written as 

\m^,my)Ai = , = — , (A19) 

2(1 + Smx,my) 



where Im^, niy) is the product of one- dimensional eigenstate of (|A2|) , and both rux and m 



"y 

are even. 

If rrix 7^ my, the Husimi function of the eigenstate is 

\{z\mx, = ^ (\{z\mx,my)\^ + \{z\my, 771^)1^ + 2Re{'mx,my\z){z\my,mx)^ . (A20) 

In the semiclassical approximation, \mx,my) and \my,mx) correspond to different tori. 
Therefore we assume 

{mx,my\z){z\my,mx) ^ 0. (A21) 
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Under this assumption, 



I {pH{\mx, my)) + pH{\my, m^))) (m^ ^ niy) 



PHi\m^,my)Aj = ^' ^'^''-^ - ' - (A22) 

PH{\mx,my)) {m^ = my) 



The second moments are 



^M2{EmjM2{E^J (m, ^ my J 



M2 = { ' ^ ^' ^ " ' ' " . (A23) 

Here, M2{E) denotes the second moment of the one-dimensional problem at energy E, whose 
explicit form is given in ( |A18D . En is the n-th eigenenergy of the one-dimensional problem, 
as given in 

If we change E^ and Ey fixing the total energy E = Ex + Ey, the product M2{Ex)M2{Ey) 
takes its minimum value at E^j. = Ey because M2{E) behaves like E" with —1/2 < a < —1/4. 
Therefore the upper limit of W2 with fixed E (we denote this as W2{E)) is 

W2iE) = (A24) 

M,(f) 

W2 of the eigenstates with m^. ~ my are close to W2, but if m^ = my, W2 of the state is one 
half of 1^2. 



2. Ergodic limit 



In this subsection, we calculate the Husimi distribution corresponding to the Berry- Voros 



hypothesis ||I2|, |13[ and the second moment thereof. 

The Berry- Voros conjecture can be stated as saying that the Wigner function for the 
stationary state of an ergodic system is approximately 



pwip,q) =NwHE-H{p,q)} 

where the normalization is, in k dimensions, 

dpdq 



w 



:5{E-Hip,q)} 



(A25) 



(A26) 



The Husimi function p^^x is obtained by smearing ( |A25|) by the Gaussian e~^^^^P'^^^. From 
the discussion in the previous section, we can see that this is equivalent to using the modified 
Hamiltonian 

k 



(A27) 
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and putting A = 1 in the Gaussian. 

We use coordinates {E, 61,62, ^3), which are related to original coordinates by the follow- 
ing equations: 



I2E 



Px = y — COS6I1COS6I2, 



Py 



2E 

— cos ^isin 62, 



where 



E^'^^ .1/2. 
^ ^ ffo \ t'lcosfc'a, 



/(^) = |i(cos^^ + sin^^) - fccos^^sin^^l 



1/4 



The normalization constant of the Wigner function Nw is determined as 

dpdq 



AttVe r/-^ d6 



I 

Jo 



V2E I ll + k 
-K 



Here, K is the complete elliptic integral. The Husimi function in this case is 

Nh 



ph{P: q) 



exp 



(A28) 
(A29) 
(A30) 
(A31) 

(A32) 

(ASS) 
(AS4) 
(ASS) 

(AS6) 



\gradH{cr) 

where Nh — N^/^fKU is the normalization constant, a is the coordinate which parameter- 
izes the equi-energy surface, and cr and ^ are orthonormal coordinates. The second moment 
of ph is 



{2Tiny 



(27r^)2V 2 



da 



86 



dO 



\gradH{e)\^ 



(AS7) 

(ASS) 

(AS9) 



Here, 



\gradH{d)\^ = 2E cos^ Oi + 

— — sin^ 6^ |cos^ ^3 + sin^ 6^ + K{K - 2) cos^ 6^ sin^ ^| , (A40) 
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and 



dcr 



dO is the volume of a small three-dimensional region formed by (d^i, d^2, f^^s) ■ After 
lengthy, but straightforward calculation we obtain 

da 



de 



Finally we obtain the following formula 



cos' e, + ^rj^ cos' e, sin' e, {f(esf + f(es)'} . 



(A41) 



M2{ph) 



where 



fZTV fix 

/ de^ / 

Jo Jq 



cos OldO 



/(^3) V2AE cos2 ei + ^giOs) sin' 9^ 



-, (A42) 



cos^ + sin^ + k{k- 2) cos^ sin' 

7W ■ 



(A43) 
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FIG. 1: Poincare surface of section at £^ = 1 and y = for = 0.0 (a), 0.2 (b), 0.4 (c) and 0.6 (d). 
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FIG. 2: Results of numerical calculation of W2 at k = 0.0. The solid line shows the ergodic limit 
calculated based on the Berry- Voros hypothesis. The broken line shows the semiclassical upper 
limit based on the torus quantization. See Appendix 




FIG. 3: W2 at A: = 0.2. 
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FIG. 4: W2 at A; = 0.4. 




(a) X=4.0 



1000 



500 



t 



++++++ 



^ + + + + + 
+ ^++++ + ^+ + + + + + - 



t + + + + + 



100 



200 



300 



E 

(b) 7^7.0 



1000 



500 



+v + +++++- 

++ * + ,+ 



+ + + + + + + + 
+ 1 ^ ^ + + + 



100 



200 



300 



1000 



500 



(c) X=10.0 



+++ ♦ ++^ + + + 

^4t-#H-+t+ + + + + + + - 



+ + 



+ + + + + + + - 
^ ^ # ^ + * + 



100 



200 



300 



FIG. 6: W2 at k = 0.2 for A = 4.0 (a), 7.0 (b) and 10.0 (c). The values of W2 depend on A, but 
the quahtative features remain unchanged. 
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FIG. 7: Comparison between quantum results and semiclassical results based on torus quantization 
at A; = 0.0. Pluses are the quantum results and circles are the semiclassical results. The broken 



line is the upper limit ( A24 ), and the solid line is exactly one half thereof. Eigenstates with 
nix = n^-y are located on the solid line. In the upper half of this figure, pluses and circles have a 
good correspondence. However, very near the broken line, the correspondence is lost because the 



assumption ( A21 ) is not good when rux ~ my. In the lower part of this figure, there are eigenstates 
with small niy, and semiclassical values are a little higher than the exact values. 
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FIG. 8: The density plot of the square of the 441st eigenfunction at k=0.2. £;=320.8, 1^2=157.2, 
NPC=102.8. The length of the sides of this figure is 15. This state is completely localized around 
the diagonal orbits, and there is no node in the transverse direction. 
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FIG. 9: NPC at k = 0.0 (a), 0.2 (b), 0.4 (c) and 0.6 (d). Dots in (b) and (c) show the regular 
series corresponding to the diagonal orbits. 
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FIG. 10: 754th eigenstate at k=OA. E=U3.6, 1^2=545.1, NPC=54.1. NPC of this state is very 
small, though W2 is not so small and its structure is not so clear. We can see scars of linear orbits 
at a; = and y = 0. 



25 



